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ABSTRACT 

The next few years will see next-generation high-power laser facilities (such as the Extreme Light Infrastructure) 
become operational, for which it is important to understand how interaction with intense laser pulses affects 
the bulk properties of a relativistic electron beam. At such high held intensities, we expect both radiation 
reaction and quantum effects to play a signihcant role in the beam dynamics. The resulting reduction in relative 
energy spread (beam cooling) at the expense of mean beam energy predicted by classical theories of radiation 
reaction depends only on the energy of the laser pulse. Quantum effects suppress this cooling, with the dynamics 
additionally sensitive to the distribution of energy within the pulse. Since chirps occur in both the production of 
high-intensity pulses (CPA) and the propagation of pulses in media, the effect of using chirps to modify the pulse 
shape has been investigated using a semi-classical extension to the Landau-Lifshitz theory. Results indicate that 
even large chirps introduce a significantly smaller change to final state predictions than going from a classical to 
quantum model for radiation reaction, the nature of which can be intuitively understood. 
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1. INTRODUCTION 

In the coming years, a new generation of high-power laser facilities (such as the Extreme Light Infrastructure 
(ELI)P|) will become operational, delivering laser intensities well beyond those currently available. The extremely 
large fields generated by these pulses will allow qualitatively new physical regimes to be explored for the first 
time, in which both radiation reaction and quantum effects will have a dominant role to play. 

An accelerating charge radiates energy (and momentum), and so must experience a recoil force. The emission 
must affect the dynamics of the particle. This is the concept of radiation reaction. Despite more than a century 
of investigation, radiation reaction remains a contentious area of physics. The most widely-accepted classical 
description of radiation reaction is based on (independent) attempts to include the effects of radiation emission by 
LorentJ^ and AbrahanP at the start of the twentieth century, which were later made fully relativistic by Dirac.l^ 
The Lorentz-Abraham-Dirac (LAD) equation for a particle of mass m and charge q in an electromagnetic field 
F reads 

-a ^ /Set ^ _±pa^^b 

m m 

where /“^t = —qF°'i,i^ is the Lorentz force. Here, the constant t := tf' j^-Km is the ‘characteristic time’ of the 
particle (~ 6 x 10“^^ s for an electron). An overdot denotes differentiation with respect to proper time. Indices 
are raised and lowered with the metric tensor r\ = diag(—1,1,1,1), and repeated indices are summed from 0 to 3. 
The i-orthogonal projection A“h := + x°'ib removes any i-component from the jerk term x, since = 0 

due to the normalisation x'^cca = — 1- This ensures that x is orthogonal to x, such that the mass shell condition, 
= —m^ where = mx^ = is preserved. Heaviside-Lorentz units are used with c = 1. 

Due to the presence of the third derivative x\ the LAD equation exhibits ‘runaway solutions’. These patho¬ 
logical solutions, which are highly unphysical and not observed, can be prevented at the cost of introducing a 
dependence on all future forces, known as preacceleration. Despite only forces within the small time r having 
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a significant influence, this response to all future forces that will act on the particle is not compatible with the 
concept of causality. The recent revie-w^^ describes these problems and proposed solutions in greater detail. 

The moset widely-used alternative to the LAD description was introduced by Landau and Lifshitz.l^ The 
radiation reaction (self-) force is treated as a perturbation about the Lorentz force, 

x“ =^ -f 0 (t), (2) 

m mm 

where once again the Lorentz force is substituted into the final term on the right-hand side. Collecting terms to 
leading order in r we arrive at the Landau-Lifshitz equation: 

. (3) 

m m \ m / 

Note that A“cF‘^^Xb = F°-^Xb + x‘^F‘^^XcXb = F°-^Xb- 

This classical equation has found widespread use amongst the community. While there is mounting evidence 
to suggest that equation @ is valid provided only that quantum effects can be ignoredpUlthe extreme conditions 
expected at ELI will take us into a regime where this is no longer possible. Instead, we must try to include some 
of the effects of quantum emission in order to make predictions for future experimental outcomes. 

Quantum effects are typically considered to be negligible if the electric field observed by the particle E is 
much less than the Sauter-Schwinger critical fieldHE3 typical of QED processes, Es = 1.32 x 10^® V/m. That is, 
when the quantum nonlinearity parameter is small. 


X := 


E 

Es 


< 1 . 


(4) 


Parameters obtainable at ELI are expected to approach y ~ 0.8, where quantum effects cannot be ignored. 
However, provided paramters are used such that ^ 1 instead, a semi-classical modification to ® should be 

valid.El 

In order to model the weakly quantum regime, we briefy consider an important difference between classical 
and quantum emission. In the classical case, the particle radiates arbitrarily small amounts of energy at all 
frequencies. However, the quantum description requires the particle to emit entire photons of energy. The 
frequency of the emitted photons is therefore limited by the energy of the particle, which introduces a cutoff to 
the radiation spectrum. The continuous classical emission which underpins the classical descriptions of radiation 
reaction is therefore expected to overestimate the radiation produced as quantum effects become important.^ 
The Landau-Lifshitz equation is therefore extended to include this reduced emission, following Kirk, Bell and 
Arka,^ by scaling the radiation reaction force by a function of the quantum nonlinearity parameter y. The full 
expression for g{x) involves a non-trivial integral over Bessel functions of the second kind. Instead, we use an 
approximation introduced by Thomas et al\^ based on fitting to numerical data, 


g (x) = (1 + 12x + 31x' + 3.7x") . 


(5) 


Predictions using this semi-classical model can be compared to those obtained from the classical Landau-Lifshitz 
theory in order to explore the role of quantum effects in the weakly quantum regime which will be available at 


ELL 


2. COLLISION WITH INTENSE PLANE-WAVE LASER PULSES 

For a linearly polarised plane wave, the electromagnetic field tensor E depends on spacetime only through the 
phase 4> = ujt — k ■ X, and takes the form 


^F\ = a{(^){e^kb-k\b), 


( 6 ) 


where k = {uj,k) is the (null) propagation direction of the laser and e the othogonal (transverse) polarisation 
vector. The function a{(j)) is a dimensionless measure of the electric field strength. An A^-cycle pulse modulated 
by a sin^-envelope has been used, for which 


o(0) 


oo sin((/)) sin^ {tt 4>/L) ioi 0 < (j) < L 
0 otherwise 


(7) 


where Og is the dimensionless (peak) intensity parameter (or normalised vector potential) and L = 2TTN is 
the (total) pulse length (full-width half-maximum duration is L/2). This pulse shape offers compact support, 
allowing the particles to begin and end in vacuum. 

The total fluence (energy per unit area) of the pulse is proportional to 

£ = I a\ct>) = ''-Nal (8) 

In this work, £ is kept constant, which fixes uq for each JV. The classical Landau-Lifshitz prediction for the 
hnal-state particle distribution is completely determined by the fluence EHEni By contrast, quantum effects are 
expected to depend directly on the intensity of the field itself, that is on the value of oq. Maintaining constant 
fluence has the advantage that quantum effects can be studied while maintaining the same classical prediction. 

Since our laser pulse is described by a plane wave, we consider the initial momenta to be strongly peaked 
around zero in the transverse directionJi^ and focus on the longitudinal properties of the distribution. The initial 
electron beam is then taken to be a Maxwellian distribution for the longitudinal momentum p (in units of me): 


f('/> = 0,p) 



(P-P)^' 

20 


(9) 


where 9 is the variance of the distribution about its mean p, and Np is the number of particles used to represent 
it. Instead of randomly sampling the distribution using a large number of particles, an iterative technique is 
used which facilitates efficient and accurate reconstruction of the distribution using far fewer particles.!^ For 
the case of classical radiation reaction according to the Landau-Lifshitz equation, this method has been shown 
to be in excellent agreement with the analytical solution of the Vlasov equationii^Hill 

In order to describe the properties of the distribution, we introduce the relative momentum spread and the 
momentum skewness, 


a{(j)) 
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respectively. The former gives a measure of the beam quality, while the latter indicates how symmetric the 
distribution is about its mean. 

We now consider the interaction of an energetic electron beam with a high-intensity laser pulse. Parameters 
have been chosen to be within the capabilities of ELL An electron beam represented by Np = 401 particles with an 
initial momentum spread of 20% around p^ = 2 x 10^ (which corresponds to an energy of approximately 

1 GeV) is collided with pulses of length N S [5,200] cycles with constant Na^ = 9248. For V = 20 with a 
wavelength of A = 800 nm, this corresponds to a full-width half-maximum pulse duration of 27 fs with peak 
intensity 2 x 10^^ W/cm^. 

Figure presents a comparison of classical and semi-classical predictions for the cooling of a relativistic 
electron beam in the interaction with two pulse lengths, N = 10 and N = 100. For a wavelength of A = 800 nm, 
the former represents a pulse of total duration 27 fs with peak intensity I = 3.96 x 10 ^^ W/cm^, while the latter 
a total duration of 270 fs and I = 3.96 x 10^° W/cm^. These intensities should be well within the capabilities 
at ELI without the need for extreme focussing. These unfocused pulses are better described by a plane wave. 

Figure [^a) shows the initial and final particle distribution for both pulse lengths predicted by both classical 
and semi-classical models. It should be noted that the classical final distributions sit directly on top of one 
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Figure 1. (Colour online) Comparison of classical and semi-classical predictions for the final-state properties of an electron 
beam colliding with pulses of = 10,100 cycles. Part (a): initial and final electron distributions. Classical predictions for 
both = 10 and 100 lie on top of one another. Parts (b)-(d): Evolution of the mean momentum p, relative momentum 
spread a, and mometum skewness during the interaction. Classical results for the final distribution show beam cooling 
with a dependence on the total fluence only, while the quantum model predicts reduced cooling dependent on the pulse 
intesnsity. 


another, supporting the statement that the final state classical predictions depend only on the total fluence of 
the laser, and not how the energy is distributed within it. The final distribution is significantly more sharply 
peaked around a much lower momentum than the initial distribution, demonstrating beam cooling. It also 
shows signs of antisymmetry, with a longer tail extending to lower momenta. In contrast, the semi-classical 
model predicts a different final distribution for each N. For N = 100, the dimensionless intensity parameter 
oo — 9.6 and we see that quantum effects have become important. Beam cooling has been reduced, and this 
trend continues to A^ = 10 (where Og — 30). 

The fluctuation of the mean momentum as the beam passes through the pulse is shown in Fig. [^b), with the 
relative momentum spread presented in Fig. [^c). For the classical cases, we see that, despite the different pulse 
lengths behaving slightly differently during the evolution due to the different numbers of cycles, both converge 
on the same final state with p = 1197.70 and a = 12.46%. For comparison, the analytical solution of the Vlasov 
equation with the Landau-Lifshitz theory for radiation reactioiP^ predicts p = 1197.62 and a — 12.53%. Instead, 
the semi-classical model for N = 100 predicts p = 1336.28 and a = 14.99%. The results for a show a reduction 
in beam cooling (d does not reduce as much from the initial 20%) with the beam remaining more energetic. 
For N = 10, the semi-classical model predicts p = 1514.40 and a = 17.36%, continuing the reduction of beam 
cooling. Note that, for N = 100, we measure the average quantum nonlinear parameter (y)^ < (y^) < 1 x 10“^. 

‘Reduction in beam cooling means that there is less contraction of phase space than would be observed in the classical 
case. As such, the relative momentum spread remains larger than that predicted by the Landau-Lifshitz theory. 





















Figure 2. (Colour online) The dimensionless intensity profile for = 10 pulses with chirp rate A = ±0.5, along with the 
original unchirped case. The pulse length La is given by equation (111. 


For the more intense pulse with N = 10, this increases to (x^) < 0.12 and we expect the semi-classical model to 
remain a valid quantum approximation. 

The momentum skewness of the distribution is given in Fig. m , which shows the initially Gaussian {S = 0) 
profile becoming increasingly negatively skewed as the evolution proceeds. The analytical solution of the classical 
Vlasov equatioiP^ predicts a final skewness of S' = —0.498, while the numerical approach used here gives S = 
—0.46^ This is reflected in Fig. [^a) by the extended tail of the classical final distribution to lower momenta. 
This is a consequence of classical beam cooling, since it is the more energetic particles which are most influenced 
by radiation reaction. The semi-classical model suppresses radiation reaction more for precisely those particles, 
and as such the distribution becomes less negatively skewed. 


3. INTERACTION WITH CHIRPED PULSES 

The use of chirped-pulse amplification (CPA) to achieve high-intensity pulses makes their influence on beam 
dynamics directly relevant to future experiments. In addition, chirps naturally arise in the propagation of laser 
pulses in media. The results of the previous section indicate that, unlike the classical Landau-Lifshitz theory of 
radiation reaction, the semi-classical model is sensitive to how the energy is distributed in the pulse, not just the 
total energy. To investigate this further, we include a chirp in our plane wave model for the laser pulse. 

We introduce the pulse length for an iV-cycle chirped pulse 


2ttN 

^ ~ 1 ± A/2’ 

where A is the chirp rate. The linearly-chirped phase is then defined to be 


li uj = d(j)/dt is our unchirped frequency, then we find 


r!(</.;A) 





( 11 ) 


( 12 ) 


(13) 


indicating that the phase contains a linear chirp. Thus, at the end of the pulse (j) = La the frequency is 
r2(L; A) = w(l ± A). (A positive chirp A > 0 therefore exhibits an increase in frequency as we move through 
the pulse.) The description of the pulse given in equation Q is then replaced by 


^This numerical difference is due to the finite number of particles used to represent the distribution.^ 
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Figure 3. (Colour online) The influence of chirping in the semi-classical model of radiation reaction. Comparison of (a) 
initial and final distributions, and (b) evolution of the spread a, with and without a chirp. The classical final distribution 
and spread are insensitive to the chirp, despite a small variation during the evolution. 
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otherwise 


(14) 


The electron beam encounters the same frequency at the front of the pulse for both positive and negative chirps, 
with the observed frequency then increasing (positive chirp) or decreasing (negative chirp) towards the rear of 
the pulse. The inclusion of the factor -^/l + A/2 = yj Lq/La in equation (14) ensures that pulses with the same 
N and oq contain the same fluence, regardless of the chirp rate. This is important, as we want to investigate 
the role of the chirp itself. A comparison of pulses of different energy, not just different distribution of energy, 
would obscure this. 


The pulse shape a{cj)]A) is shown in Fig. for A = ±0.5, along with the unchirped pulse (-). This 

illustrates how the increase in frequency caused by the positive chirp ( - ) generates a shorter pulse length 

than the unchirped pulse. The amplitude is increased from oq (indicated by the horizontal dotted lines) such 
that the fluence remains constant. The negative chirp (— ■ — ), on the other hand, causes an extended pulse of 
lower peak intensity for the same fluence. 


Since the fluence is kept constant, predictions from the classical Landau-Lifshitz theory for the final properties 
of the beam are not affected by a chirp. This is confirmed in Fig. where the influence of a chirp in the semi- 
classical model is also shown for the interaction of a 1 GeV electron beam with N = 20 (ag — 22) laser pulses. 
Figure l^a) shows the initial and final distributions obtained using both the classical and semi-classical models, 
with and without chirps with A = ±0.5. The evolution of the spread a as the bunch travels through the pulse 
is shown in Fig. [^b), which highlights the significant reduction in beam cooling when using the semi-classical 
model. Even with a relatively large chirp of |A| = 0.5, the impact of chirping the laser pulse can been seen to 
be a notably smaller effect than moving from the classical to the semi-classical model. 

For the positively chirped case, a final spread <7+ = 16.9% is found, compared to do = 16.6% for the unchirped 
pulse and d-_ = 16.3% for the negative chirp. These changes in the spread are easily understood by considering 
the peak intensity encountered in each case. For the positive chirp, the pulse length is reduced because of the 
decrease in wavelength as we move through the pulse. To keep the fluence the same as for the unchirped pulse, 
the peak intensity must increase from ag to ao\/l ± A/2. In turn, this increase leads to a higher instantaneous 
value of the quantum parameter, y, and thus increased suppression of the radiation reaction effect. Hence, beam 
cooling is further reduced by a positive chirp. In contrast, a negatively chirped pulse has a longer duration and a 
corresponding lower peak intensity, therefore quantum effects are less important and the radiation reaction term 
is less strongly suppressed. 
































4. CONCLUSIONS 


In the next few years, as new high-powered laser facilities come online, it is increasingly important to understand 
the physics that will occur at these unprecedented field strengths. Previously untested areas of physics will be 
experimentally probed for the first time, which requires a knowledge of fundamental principles such as radiation 
reaction in regimes for which quantum effects can no longer be ignored. 

In this paper, a semi-classical extension to the Landau-Lifshitz theory of radiation reaction has been used 
to investigate the role of chirps in the interaction of relativistic electron beams with high-intensity laser pulses. 
Comparison is hrst made to classical predictions to identify the differences as some quantum effects are included 
in the model. Classical results are compared to analytical solution of the Vlasov equation, supporting the validity 
of the numerical technique used. A reduction in the amount of beam cooling is observed in the quantum case 
and found to depend on the intensity profile of the pulse itself, unlike classical predictions which only care about 
the total energy of the pulse. Chirps have been used to alter the distribution of energy within the pulses, with 
even a large chirp of 50% resulting in only a small modihcation to the hnal state properties, the origin of which 
is intuitively understood. 

The results presented here are limited to the semi-classical case <C 1 due to the deterministic equation of 
motion used to describe radiation reaction. Strongly quantum regimes may be explored using similar techniques 
by adopting a stochastic equation where photon emission probabilities are determined by strong field QED, as 
inEHEH This will be addressed in future work. 
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